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SUMMARY 

SmaLL-perturbation spectral-analysis techniques are used to obtain 
the root -mean- square flame- generated turbulence velocities and the at- 
tenuating pressure fluctuations stemming from interaction of a constant- 
pressure flame front with a field of isotropic turbulence in the absence 
of turbulence decay processes. 

The anisotropic flame- generated turbulence velocities are found to 
be of about the same intensity as those of the incident isotropic turbu- 
lence, the lateral turbulence velocities being always lower but the lon- 
gitudinal velocity is somewhat increased for flame- temperature ratios in 
excess of 7. The small-perturbation analysis indicates that the incre- 
mental turbulent flame speed is a second-order quantity con^osed of two 
parts: One part represents the root-mean- square area of the turbulent 

flame front; the other represents the contribution of the transverse ve- 
locity fluctuations which resiilts from the flame-front distortion. Di- 
rectly at the flame front the noise -pres sure levels of the pressure 
fluctuations are fairlj intense (59 to 81 decibels referred to 0.0002 
dyne/sq cm) even at moderate approach- flow turbulence intensities. 


INTRODUCTION 

Development of hi^-output Jet engines has stimulated interest in 
the role played by turbulence in combustion phenomena. In the earliest 
studies of flame- turb\olence interaction, Damkohler (ref. l) and Shelkin 
(ref. 2) utilized mixing-length theories of turbulence to obtain semi- 
quantitative relations for predicting flame speeds. Damkohler intro- 
duced the concept that turbulence of a scslLe large relative to the 
flame-front thickness increases the average flame speed by increasing 
the instantaneous flame surface area. The authors of reference 3 were 
unable to confirm the relations of references 1 or 2 by their experi- 
mental values of turbulent flame speed in Bunsen burners. Experiments 
on flames stabilized in channels (ref. 4) suggested that approach- flow- 
turbulence iiad little effect on burning velocity and that the disturb- 
ances affecting turbulent flame speed were primarily flame-generated. 

A similar conclusion was arrived at in reference 5. 
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In an attempt to obtain agreement between theory eind experiment, 
Karlovltz, Denniston, and Wells (ref. 6) arid Scurlock and Grover (ref. 

7 ) have incorporated the concept of flame^enerated disturbances in 
their recent theories of turbulent flame speed which utilize G. I. 
Taylor's one -dimensional theory of diffusion by continuous movements. 

The somewhat arbitrary assumption is made in the analyses of both ref- 
erences 6 and 7 that these flame -generated disturbances constitute ad- 
ditional turbulence only. In reference 6 the energy of the flame- 
generated turbulence is teiken as the difference between kinetic energy 
of the burned gas in the absence of turbulence and the kinetic energy 
obtained by using the average velocity of the burned gas normal to the 
turbulent flame front. In reference 7 the flame -generated turbiilence 
energy is obtained from a mcmentiom balance- of unburned and burned gaS&s 
before and after an assumed mixing of the burned gas. 

The data obtained in reference 8 on pentane-air flames baffle- 
stabilized in a rectangular duct suggest that the methods of references 
6 and 7 overestimate considerably the turhulence generated by flame - 
turbulence Interaction. Apart from the question of validity of such 
methods of calculating flame- generated turbulence, objections have been 
raised (ref. 9) to calculations of flame'^eed made on the basis of a 
hypothetical upstream turbulence compounded from approach- stream turbu- 
lence and turbiilence generated downstream' of the flame as was done in 
reference 7 and implied in reference 6. 

The present analysis is primarily concerned with the turbulence ve- 
locities and other fluctuation quantities associated with the linearized 
interaction of a free flame (not influenced by bounding walls) with tur- 
bulence present in the combustible mixture. Such turbulence will be re- 
ferred to as approach-flow turbulence. The flame is treated as a dis- 
continuity specified by the appropriate fundamental (laminar) flame 
speed and flame temperature. The interaction of such a flame front with 
a transverse plane wave, that is, a vorticity wave or shear wave, of ar- 
bitrary inclination relative to the front is first analyzed. The ef- 
fects of an entire . spectrum of transverse plane waves constituting a weak 
field of turbulence are then developed from the* single-wave results. The 
statistical or root -me an- square fluctuation qtiantities describing the 
pressure fields and the anisotropic flame -generated turbulence resulting 
from interaction of the flame front with isotropic approach- flow turbu- 
lence are obtained for the limiting case of constant-pressure combustion. 
Some discussions of tiirbulent flame speed and of combustion noise are 
also presented. 


FLAME - TURBULENCE INTERACTION PROCESS 

Turbulent motion may be regarded as a Foiirier superposition of a 
very large number of different- si zed and randomly oriented component 
plane-wave motions. The customary assumptions concerning the turbulence 
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(see ref. 10) are also made in the present analysis - namely, that tmr- 
■* bulent decay effects are ne^igible and that the density fluctuations 

associated with the turbulence are also small enough to be neglected. 

The first assiimption, which iit 5 )lies inviscid flow and very small turbu- 
lent velocity fluctuations, permits linear superposition of the compo- 
nent waves. With the second assumption, the continuity equation re- 
quires that these Fourier waves be transverse plane waves, that Is, 
vorticity or shear waves. For each of these waves the local velocity 
vector A is perpendicular to the vector k, normal to the wave front. 
The vector k is termed the wave-numdier vector; its magnitude k is 
termed the wave number, which is defined as 2it divided by the wave 
length. All symbols are defined in appendix A. Any one of the parallel 
planes containing both the local velocity vector A and the wave-number 
vector k is called the "polarization plane." 

Because of the assimned linear superposition of the conponent waves, 
the ccmplete Interaction results can be obtained from study of the in- 
teraction of a plane flame front with a single- component transverse wave 
of the turbulence field. For simplicity this typical vorticity wave 
<1 will first be taken as a two-dimensional wave. Generalization to the 

three-dimensional case will be made later. The configuration considered 
is shown in sketch (a); 
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A flame front moves 'witb. mean velocity U|p inta an inviscid combustible 
mixture. This mixture Is at rest, but contains a vortlclty wave with 
velocity vector A (components and Vg and with wave-nimiber 

vector k Inclined at an angle <p to the positive direction of the 
Xg-axls. In the absence of any perturbatlcm Interaction, the plane 
front propagates Into the combustible mixture with velocity U (the 
laminar flame speed). As Indicated In appendix B, the motion of the 
bxxrned gas whose velocity V Is constant Is away from the flame front. 
The flame front Is assumed to be of Infinite extent In directions trans- 
verse to the direction of the x^-axls. 

References 11 and 12 have shown that a weak Invlscld disturbance 
field may be resolved into a stationary component and a moving cotDpo- 
nent, both relative to the mean local flow. The moving component is 
an irrotational isentropic pressure-velocity disturbance. The station- 
ary component, which is convected by the mean, local flow, is a constant- 
pressure disturbance containing any vortlclty fluctuations and entropy 
fluctuations present in the disturbance field. Thus, the interaction 
of the flame front with the vortlclty wave would be expected to generate 
both an irrotational disturbance with velocity components ^ ^p b 

and a rotational disturbance with velocity components Ug -jj, Vg ^ In 

the burned gas, and an irrotational disturbance with velocity components 
Up^c^ ■''■p c combustible mixture. The resulting velocity fluctua- 

tions, which include both the irrotational and rotational disturbances, 
are designated as u^, v^ and u^, V'^ for the ccmbustible mixture and 
the burned gas, respectively. The flame front is displaced by an amount 
C(y^'t) from its mean or unperturbed position as a result of the 
interaction. 

The diagrams of figure 1, which are si^lar to those used in refer- 
ence 12, may prove helpful in visualizing the interaction process. Sup- 
pose that at some , instant t^_ the flame intersects a front of the vor- 
ticity wave at point of figure 1(a). At a later time t]_ + Bt, the 

flame has moved a distance UBt and now intersects the stationary vor- 
ticity wave front at point P 2 . A vorticity wave with front parallel to 
line QPg is then produced in the b-urned gas. A cylindrical sound wave 
is generated at point P 2 _ at time t^ and propagates at speed a^ 
into the burned gas while being convected with velocity -V. Another 
cylindrical wave is generated at point Pj at time t^^ and propagates 
into the combustible mixture with speed a^. The cylindrical wave 
fronts thus generated form envelopes (Mach lines) in both the combusti- 
ble mixture and the burned gas which constitute plane sound waves. — 

For the wave- inclination angle cp shown in figure l(b), an enve- 
lope cannot be formed on the burned-gas side of the flame front. The 
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cylindrical sound waves thus expand independently and are thereby atten- 
uated. On the combustible-mixture side of the front, the cylindrical 
sound waves meet at the common tangent point Pg. For inclination an- 
gles less than the critical angle shown, attenuating pressure waves are 
also produced in the combustible mixture’ until another critical angle 
180° - <P is reached. Below this second critical angle, plane soiind 
waves are again obtained. These critical angles may be obtained from 

the geometry of figure l(b). For small flame Mach numbers where M= 

^c 

‘Pcr,c = sin-1 M 




^cr,b 



As the flame Mach number M decreases, attenuating waves are produced 
for a wider range of inclination angle. In the limiting case of very 
slow flow (constant-pressure combustion),! only attenuating pressure 
waves appear in combination with the vorticity waves if 0° < (p < 180°. 
Quantities associated with the pressure wave vanish when cp = 0° or 
180° because the Incident vorticity wave then passes throu^ the com- 
bustion front without distorting the front. 


SIlfGLE-WAVE ANALYSIS 
Two-Dimensional Formulation 

The interaction process described in the preceding section is now 
formulated analytically for the passage of a combustion front through 
a single weak two-dimensional vorticity wave of constant density in- 
clined to the flame front. The case of a vorticity wave in three di- 
mensions is considered later. 

The combustion front is assumed to be ccm^^letely specified by its 
laminar flame velocity U and the ratio of stagnation temperatures t 
in the burned gas and in the combustible mixture, respectively. In 
the absence of any perturbations, the equations for conservation of 
momentum, energy, and mass-flow rate, respectively, as written for a 
reference frame moving with the flame front, are 


■^It can be shown that the static-pressure ratio across a flame 
front is given by P-b/p,, = 1 - - 1 )m2 + . . . 
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Pc + = Pb + Pb^U - 

( 7 ^ ''c + I U 2 ) = ^ Tb + I (U - V)2 j> 
PcU = p^,(U - V) 


( 1 ) 


vhere subscripts c and b designate stations in the combustible mix- 
ture and in the biorned gas^ respectively. For sinplicity it has been 
assumed that the combustion process does not entail a change in the 
number of mbles per unit mass of gasj also, differences in the ratio 
of specific heats for the burned and unburned gas are ignored. The 
quantity ('t - l) is then indicative of the increase in stagnation en- 
thalpy or heat release. 


For the interaction problem the resulting flame-front distortion 
S'Cyj’t) must be considered in addition to the generated disturbemces 
previously mentioned. Thus, both the flame-front perturbation veloc- 

ity ■ 3 ^" ^ ^ PtL® Instantaneous flame-front slope - Ky vill 


appear in the equations of motion. The conservation equations may still 
be applied in a coordinate system moving instantaneously with the dis- 
torted flame front since extreme gradients occur across the front and 
small disturbances are postulated. The various perturbation quantities 
(designated by primes) are assumed to have zero space or time averages. 
The flame speed Uij (see sketch (b)) will thus Include any time- 


independent contributions arising from the perturbations. 
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Conservation of normal and tangential momentum, energy, and mass-flow- 
rate provide -fclie following relations: 


(Pc + + tPc Pk) [(“r + n - ^k^y^ * 2(1^ + Vi - u;)vie'] 

- (Pb + PbHl + + (Pb Pb) [("t - - 1^)2 + 2(u^ _ T,. + t; - 

+ £t “ ^c)^y “ ‘'^c " * Vgj + £t “ ^ 

(t + (To * (% + 5t - ■*'c^ ■"7^ (Tb + T^) + I (1% - + €t - “b)^ * I 

(Pj. + Pq)(^ ■•■ - “c ■*■ " “b ■•■ 


(2a) 

(Zb) 

(2o) 

(2d) 


Small-perturtation techniques are used to make the interaction 
problem amenable to analysis. Haen, if the velocity perturbations are 
assumed to be small relative to the flame speed % and the flame-front 
slope is also assumed to be very Rjna .1 3 so that terms of second or- 

der in the fluctuation quantities may be neglected, application of the 
linearized state equation = ~~ + utilization of equations 

(l) permit the following boimdaiy conditions at the flame front to be 
obtained from equations (2): 


Pb 

Pb 

Pb ^-t " ^b -r, „ ^t " “c 

+ Bi ~ + B 2 u B 3 + B 2 u . 

(3a) 


+ Ei^y 

(3b) 

Pb 
Pb " 

Pb ^t - ^b ^ ^ - ’^c 

(3c) 


p; ^ ^i-% ^i-K iK 
u u r Pc 

(3d) 

- _ ru^ fu - V 
a|U 

T- >»■ 4) 


1 

III 

(^), . 11 ^. 1 K,. - 1 ) (r . ^ 

a| = ajj = ^ 

^ (4) 
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It has also been assxuned that the flow upstream of the flame is isen- 
tropic. It will be shown In the section Turbulent Flame Speed that 
Uy = U is correct throu^ first-order terms. 

Another relation is required at the flame front. For the two- 
dimensional case under consideration, the local instantaneous normal 
propagation velocity U + 5U of the distorted flame front into the 
combustible mixture at rest (see sketch (b)) is 


U + 8U = 


U, 


T 


£ ‘ - 


VL, 


r'X 


4 


1 + 


£•2 

y 


The incremental propagation velocity 5U will be determined from exist- 
ing Information on laminar flames. Seme support for this procedure is 
given in reference 13. As reported therein, radiant flux- intensity 
measurements on laminar and turbulent prqpane-alr flames suggest that 
a small surface element of a turbulent flame is chemically and physi- 
cally the same as that of a corresponding laminar flame. 


The propagation speed of a l amina r flame is affected by both the 
ambient pressure and the ambient temperature. Althou^ the functional 
relations have not been rigorously determined, preliminary indications 
are that the pressure effect is much smaller than the tenrperature ef- 
fect. In the present anetlysis the locad. flame propagation speed is 
asstmed to depend only upon the ten^jerature of the combustible mixture, 
that is, 


8U 




With the en 5 )irical relation obtained in reference 14 as a guide, it is 
assumed that U = r^ + where r^, v^, and n are constants which 

depend upon the fuel and oxidant under consideration. Thus, 




and 
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Correct throu^ first-order terms -wherein Hj> = the following bound- 
ary condition is obtained at the flame front: 


^t 


u. 


U 



(3e) 


In a coordinate system fixed in space, the equations of motion for 
the -two-dimensional fluctuation quantities in the burned gas -with terms 
of second order neglected, are 



Sui 



J_^b 



'ST 

+ V 


' Pb ^ 



^b 


^b 




"ST 

+ V 

bx^ 

Pb 




Spb 



1 n 


+ V 




1 = ° 





f^b 


*^b°V St’ 

+ Pfc 

)®V^ 

dxc 


dy J ~ 


For a coordinate system moving -with constant velocity V, the preceding 
equations reduce to the same form as the corresponding equations for the 
fluctuation quantities in the combustible mixture relative to a coordi- 
nate system fixed in space. Thus, the flow equations for both the 
combustible -mixture fluctuations and the burned- gas fluctuations may be 
written, with appropriate subscripts c or b, as 

Bu^ _ 1 I 

^ " p ^ 


Sv ' _ 1 3p ' 

dt ~ p Sy 



ST' 

'ST 




( 5 ) 
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The coordinate Bystems y and x-jj, y for the first-order analysis 

are indicated in sketch (c). 

t=0 



A two-dimensional vorticity wave in the combustible mixture with 
velocity vector of magnitude A and having its wave-number vector k 
inclined at an angle <p to the positive direction of the X(,-axis may 
be written in the form 


u 


U 


= (A sin <p)e 


iv 


s,c 

U 


= (- A COB <p)e^^ 


where v s k 2 _Xg + kgy, and k^[ and kg are camponents of the wave- 
number vector k in the Xg- and the y-directions, respectively, with 
kg/k^ = tan <P. As a result of the linear boundary conditions of equa- 
tions (3), obtaining a imlque solution of the interaction problem re- 
quires that the arguments of all disturbance waves match at the flame 
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front. This matcMng requirement together with provision for differ- 
ences in phase angle yields the following fomi for the vorticity wave 
present in the bnmed gas: 



(Gi + 



+ il2)e^+ 


where 







Pressure fluctuations generated by the Interaction must satisfy the 
following wave equation, with appropriate subscripts c or b, which is 
obtainable from equations (5): 





= 0 


( 6 ) 


The present analysis will be concerned with the limiting case of very 
slow flow (constant-pressure combustion). It is clear frcm the rela- 
tions given for the critical wave-inclination angles <Pcr,c ‘Pcr,b 

that, for very slow flows, only the irrotational isentropic pressure 
waves described in the section FLAME - TURBULENCE INTERACTIQR PROCESS 
will be generated by the Interaction. The form of these pressure waves 
that satisfies equation (6) has already been established in reference 12 
in terms of the variables t] and 5, where 

Tic - (^‘c)(*C - 

?c - ^c^c + c^y + ua^,t 

r ( 7 ) 

lb = 

The variable ti Is prcportional to the distance from the flame front. 

At the front, = Oj upstream of the flame, is positive; and 

downstream of the flame, t}^ is positive. The equation ^ = constant 
defines planes moving with constant velocity {Ud)^ at an angle 
tan“^(c/b)j, to the flame velocity U. Equation (6) takes the form 
of Laplace's equation 
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In ^ = «2 |n.|I 



( 8 ) 


Matching arguments of the pressure and vortlcity vaves at the combustion 
front vhere Tj = 0 and satisfying the requirements of equations (8) 
provide the following values for the constants of equation (7): 


c 


c 




<ib = 








( 9 ) 

In addition to the boundary condition from, equations (3), the pressure 
fluctuations will be required to satisfy the boundary condition — = 0 
at T\ = ». Utilizing equations (36) and (37) of reference 12 yields 




= (R^ + lR2)e^^c-’lc 

(10a) 

= (Ji + iJ2)e'’^b"''b 

(10b) 


In the combustible mixture, density fluctuations are associated only 
with the presstire fluctuations according to the isentropic relation 

Pc 1 Pc 

In the burned gas, density fluctuations may also be caused 

*^c * Pc 

by entropy fluctuations generated by flame distortion and heat-release 
fluctuations, if present, as well as by pressure fluctuations. Velocity 
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fluctuations are asaociated with both pressure fluctuations and vortic- 
ity fluctuations. It is convenient to deal with the pressure coeffi- 
P' 


cient 


P.U2' 


Thus, the disturbances arising from the interaction of 


the flame front and incident vorticity wave take the following forms: 

_ I 


PcU" 


II 

(Hi + 

lR2)e^^c-'lc = 

(r(1) 


(10c) 

-A 

' PcU' 

3 (Jl • 


* (j(l) 

+ ij(2))e^^^~’^^ 

(lOd) 

II 

-^Itf 

f (Jl 

+ iJ2)e^^l3"’^t. 

+ (L-l ■ 

+ lL2)eit 

(lOe) 

U 

- (Ni H 

- ilT2)e^^^‘’^^ - 

f (G]^ + 

i02)ei^ 

(lOf) 

u 

= (Pi -f 

- iP2)e^^^~’^^ - 

^ di + 

il2)eit 

(lOg) 

^c 

u 

= (Wi 


+ (a sin cp)e^'^ 

(lOh) 

1 

^c 

u 

= (Xi 


- (a cob <p)e^'' 

(101) 


The flame displacement velocity may be written 

^ i(kTUt+kpy) , , 

■rf- = (Hi + iH2)e = (H^ + iH2)e 


(lOj) 


To satisfy the requirements that the arguments of all disturbance quan- 

tities match at the flame front and that the flame- front 

slope must take the form 


Cy = (Hi + iH2)(tan 


(10k) 



14 


KACA TN 3407 


Integration of equation (lOj) witli respect to time t gives the flame 
displacement as 


K' (Hg - i%)ei‘^ = (Hg - iHi)ei<^ (lOl) 

J£l ^2 


Extension to Three-Dimensional Disturbances 

Equations (lO) describe the interaction of a flame front with a 
constant-density vorticity wave having velocity components 


u, 


U 


= (A sin <p)e^'' and 


■BjC 

U 


= (- A cos q>)e^'' 


in the and y-directions, respectively. The vorticity wave may be 

Wg 

considered to have a third velocity component — = Ce^'' in the z- 

direction. In the preceding linearized analysis, the amplitude C was 
not prescribed. This component, which is normal to both the u’ and 
v' components and parallel to the plane of the flame front, then is 
associated with a corresponding congjonent of the vorticity wave in the 

V * -L 

burned gas — 


Inasmuch as turbvilence fields are three-dimensional, the interac- 
tion equations must be revised accordingly for application to the spec- 
tral analysis which follows. Assume, as shown in sketch (d), that the 
polarization plane which will contain the wave-number vector k» is in- 
clined at some angle 9 to the x-j_ j,, Xg ^-plane of a new coordinate 
system Xg^^,, Xj^^, fixed in space for the combustible mixture. 




f “ 
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The correBponding coordinate system ^3,b burned 

gas is assumed to be moving with mean.- flow velocity V. Components 
ki, kg^ nnri kg of the wave-number vector k in the directions of 

^l,c^ xg^c^ ^3,c^ respectively, are 


k]^ = k cos <p 
kg = k sin (p cos 6 
kg = k sin <p Bln ej 


( 11 ) 


With primed perturbation vector quantities referring to the original co- 
ordinates Xq, x-jj, y, z and unprlmed perturbation vector quantities re- 
ferring to the xi^g, xg^c^ ^3,c ^l,b^ ^2,b^ ^3,b coordinate sys- 
tems, the following transfoimatlon relations apply: 


l,c 

= V* cos 0 - w' sin 0 


Un 


l,b - '^b 

Ug^-jj = v^ cos e - sin 0 
= Vg sin 0 + Wg cos 0 ’^3,b “ '''b ® + "^b ® 

C = €• 


^2,c ” ''c ~ "c 


^3,c 


} ( 12 ) 


«t =«t 


- «: 

^* 3,0 ■ 


COS 0 


sin 0 


J 


This notation refers only to equation (l2). 


Wi-th the use of eqxiations (lO) and (12), the interaction fluctua- 
tion quantities (again designated by primes) referred to the coordinate 
^l,c^ ^2,c^ ^3,c ^l,b^ ^2,b^ ^3,b W written as 
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PcU" 


= (r( 1) + iR(2))e^^c"^c 




P^U2 


^ = (jd) + ij(2))e^^^-^t 


— = i 

Pb 


^ i (El + lE2)e"^<='’'= 


. ^ (Jl + lJ2)e*‘^*>‘'''> + (Ii + lL2)el* 


u, 


2>e 

U 


= (Wi + iW 2 )e^^c-^c + (A sin <p)eiv 
= (X^ + iX 2 )(cos 0)6^^*^ - (a cos <P cob 0 + C sin 0)e^'' 


- ^ •~ = (X^ + iX 2 )(sin 0)e^^° - (A cos 9 sin 0 - C cos 0)e^'’ 




u. 


u, 


2,b 

U 


^,b 

U 


^ = (Hi + lN2)e^^^”’^^ + (Gl + iG2)e^^ 


(Pi + lP 2 )(cos + [(li + ll 2 )cos 0 - C sin 

(Pl + iP 2 )(Bln 0)e + |(li + il 2 )Bln 0 + C cos 

^ = (Hi + iH2)ei^ 

^X 2 2 (Hi + iH 2 )(tan <p cos 0)e^P 
= (Hi + lH 2 )(tan <p sin 

3, c 


_ (H2 - iHi)e 

^ ~ k' 


ia 


J 

(13) 
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!The fluctuation £tn 5 >lltude coefficients of equations (lO) or (13) 
may be determined from eqiiations (3) and (5). Details of the solution 
are given in appendixes B and C. The general solution for the attenu- 
ating pressure- wave regime is indicated in appendix B. Inasmuch as the 

flame Mach number M s JL is generally much less than 0.01, the limit - 

®"C 

ing case of very slow flow (constant-pressure combustion) provides a 
reasonable sin^ilifl cation of the problem and only the attenuating-wave 
solution need be considered. The amplitude coefficients for this lim- 
iting case are given in appendix C (eqs. (C8)). v 

Equations (13) and (C8) describe the linearized interaction of a 
constant-pressure flame front with a single vorticity wave or shear 
wave having its wave-number vector k inclined at an angle (p to the 
direction of travel of the undisturbed flame front and having its plane 
of polarization Inclined at an angle 9 to the x^^ j,, 

the coordinate axes ^2,c^ ^3,c* wave-nu^er vector of the 

shear wave generated in the burned gas makes an angle 

<Pb tan“l('E’ tan cp) with the direction of propagation of the undisturbed 
flame front. Atteniiating potential fields are generated in both the 
combustible mixture and the burned gas. Physical quantities associated 
with these fields attenuate e3q>onentially with increasing distance frcm 
the flame front. The amplitude coefficients for a given heat release 
(a prescribed t), a given inclination angle <p and polarization angle 
0 depend upon both the Intensity of the incident vorticity wave and the 
heat-release perturbation parameter t'/t. In the absence of such heat- 
release perturbations, there are no density fluctiiations in the burned 
gas (correct to order M^). 

These single-wave results may be used to determine the Interaction 
of a constant-pressure combustion front with a turbulence field of con- 
stant density for the case of negligible turbulence decay. The turbu- 
lence field will contain an infinite number of transverse plane waves 
with all wave lengths and planes of polarization. The spectral analysis 
technique used in obtaining such a superposition of waves will be dis- 
cussed briefly before proceeding with the interaction problem. 


SPECTEAL AMEISIS, 

Gteneral Considerations 

A turbulence field satisfying the incompressible-flow continuity 
equation may be represented., by the following superposition of plane 
transverse waves: 
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u(x,t) =yyy'e^^*2^(k,t) 

where x is a position vector, k is a wave-numher vector, t is the 
time, and ^(k,t) is the random amplitude vector of a component wavej 
and the subscript i takes on the values 1, 2, and 3. The quantity 
^(k,t)e^^*2E represents the contribution to the velocity field from a 
volume element dk in wave-number space. When, as in the present case, 
the equations of motion are linear there is no modulation or interfer- 
ence between con^jonent waves, and the various statistical quantities 
describing a randcan field may be obtained from the results of a single- 
wave analysis. To avoid the interpretative difficulties associated with 
the random variable Z(k,t), which is nondifferentiable with respect to 
k, use is made of the techniques of references 10 and 15 which utilize 
correlation spectra rather than an^ilitude spectra in the analysis of 
homogeneous tiirbulence. 

■ A velocity correlation is defined as the ensemble average 
Ui(x,t}uj(x' , t) of the product of a fluctiiation-velocity component u^ 
at X and a component uj at x' = x ■*■ £ where r is a separation 
vector. The ensemble average, designated by a bar, may be regarded as 
the result of averaging the product ■u^(x,t)uj (x,t) at a given instant 
over a very large nu mber of statistically similar fields. The nine ve- 
locity correlations constitute the velocity correlation tensor 

Tij(x,x‘,t). For a homogeneous field, T^j^j depends only on r, so that 
the tensor may be written 

As shown in references 10 and 15, the velocity correlation tensor 
has the following Fourier integral expansion: 

- (14) 

where is the spectral tensor function, ^j^j(k,t) is the spec- 

tral. tensor density of a homogeneous turbulence field and 

#^j(k,t)dk = dZ*(k,t)dZj(k,t) (15) 

where dZ*(k,t) denotes the complex conjugate of dZ^(k,t). For r = 0 
and i = J, equation (l4) may be written 
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Tii(0,t) 



dzJ(k,t)dZi(k,t) 


(16) 




CD 

W 


O 

t 

-3 


For homogeneous turbulence fields, wherein ensemble averages and ^ace 
averages are identical, equations (l5) and (16) provide the basis for 
obtaining the spatial mean-sgtiare velocity ceai^jonents from the single- 
wave results given by equations (13) and (C8). Equations (15) and (16) 
are also applicable to scalar fields. In the absence of viscosity, as 
postulated, the shear-velocity fields present in the combustible mixture 
and in the burned gas are homogeneous and application of equation (16) 
presents no complications. The corresponding potential- flow fields, 
althou^ spatially inhomogeneous, are homogeneous in the given 
X 2 ,X 3 -planes. It has been shown in reference 16 that equation (16), 
in effect, may be applied for such fields to obtain the mean- square 
fluctuations pertaining to a given plane of homogeneity. 


As a result of the preceding, discussion, the single-wave interac- 
tion results for constant-pressiire combustion will be used to obtain 
the spectral densities of the fluctuation quantities at the flame front 
where the attenuaticjn factors and e~^b are unity and 

= 0 = t = V. For conciseness define RA sin <p = r(1) + 1r(2), 
WA sin <p s Wi + iWg, GA sin <p s Gi + iGg, HA sin q> 3 Hi + iHa, 

HA sin (p = Hi + iH 2 ^ and C^A sin <p sC. With the notation 


u? = rrTdz* dz^ 


and, for example. 



as in equation (16), the following equations are obtained by analogy 
with equations (13) and (C8) for the case where heat-release pertur- 
bations are absent, ttiat is, t’ = 0: 
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= R dZ 


#) ■ #) ■ 

- IW cos S - (cot q> cos e + C^sln 0)aZj^^c 

T 

= - iW sin e dZjL^c " <p sin 0 - C^cos 0)dZ3_^c 

d(%) = iH cos 0 - (^ot ^ ^ e)dZ^^, 

d(:ii) = iH Sin 0 dZi^, - {S2L±^ G - C^cos 0)dZi,, 


} (17) 


d(«') = - ^ <^l,c - “ k cos (P “^l^c 

d(^) - H dZi,„ 

H«i2,c) = H tan tp COB 0 
d(5^^ ) = H tan <p sin 0 




In the velocity ratios of equations (17 ) , the first term on the right 

^i,p,c 

side represents the potential- flow contributions ^ or 


U 


Ui g 2 

the second term represents the shear-flow contributions y ^ 
subscript i takes on the values 1, 2, and 3. 


5493 
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At a given instant the spatial mean-square potential- shear- 
flow contributions of i the distirrbance fields are obtained from.' 
equations (17): 




-fff 










> ( 18 ) 






The taean-square flame front quantities are 




(h*h) 

,2 2 -dZ, - 

k^cos^(p 


d +€L 
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From equations (C8) of appendix C, with. ,t' = 0: 


R*R = 


(t - 1)^(t tan^q> - l)‘ 
A sin^q> 


W*w = - l)^(T^tan^(i> - 1)^ 

Q*Q _ T^sec^<p[A - 4 t(t:^ - l)tan^<^ V 
A(l + T^tan^cp) 


N*N = 


t^(t - 1)^(t tan^q) - l)^sec^(p 


H*H = 


A(l + T^tan^q)) 
2 

sec q> 


4t2.-2. 


J 


(19) 


Mean-Square Fluctuations for Initial Isotropic Turbulence 

For a given combustion process (t and U prescribed) , the spatial, 
mean- square fluctuations of equations (18) depend upon the quantity 

'^^l,c‘^l,c^ whiGh is specified by the type of turbulence present in the 

combustible mixtxire. The res-ults obtained in reference 16 for the in- 
teraction of axisymmetric turbulence with a shock wave suggest that for 
the present problem the degree of smisotropy of the incident turbiilence 
field may not be critical. Fxir simplicity, the turbulence in the com- 
bustible mixture is assumed Isotropic. 

As indicated in reference 10, the spectral density tensors for any 
isotropic turbulence field satisfying the incompressible-flcw continuity 
equation are 

$ij(k) =a(k)(k^6ij - k^kj) (20) 

where k^ = kf + k^ + q ^ for 1 = jj and 5. , = 0 for ± f y, 

-L ^ O Xj Xj 

and is the scaleir amplitude function defining the spectral density 

tensor. From equations (15), (20), and (ll). 


dZ^^i,dZ]_^Q = 2(k)k^sln^<p dk 
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or, transforming to spherical polar coordinates k, B, wherein 
dk = dkj^dkgdkg = k^sin q> dk d0 d«p, yields 


^, 0 *^ 1 , c = de sin^q> d<p 




dZ* dZp = 2(k)k'^dk d0 sin ^(cos^(p + sin^«^ sin^sldcp 

dZ* dZ<, = 2(k)k'^dk d0 sin ipfcos^p + sin^jp coB^oJdtp 

O ^ C o ^ C 


\ ( 21 ) 


dZ* dZ„ + dZ* dZ_ = 2(k)k~dk dB sin <p[l + cos^c^dcpx 

O^OO^C * 

The mean- square velocity components of the incident isotropic turbulence 
field are then given by 


— « 2jt It e» 

f- a(k)k^dk dB sin^q> dq> = ~ S(k)k^dk 

'-'q ■^0 




^ 2(k)k^dk dB JP sin (P(l 

m 


+ COB 


2q>)dq> 


> 


J 

(22) 


The remaining spatial mean-square fluctuation quantities, referred 
to the intensities of equations (22) in order that the scale need not 
be specified, are given by 


(^)' 


3 f 

4 


4K 

- l)^ - (t tan^ip - l)^dtp 


(23a) 
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gy)‘ f 

Pif? Pf4 ‘ ' “ 


(T tan2(p 


( 3t^ Fa - 4t(t^ - l)tan^<p]sln q> tan^q) 

— ^ ~ 4: J A(1 + T^tan^ep) 




d<p 


(~^^) ^ 1 _ 3 (^2 _ 1) r [a + 4t sec^qJslaSq 

— ” 8 ' J aQl + T2tan2<pJ 


(p:^ ()2^) » (-5^) _ s,2,, ^ f (t t«.^ - p;ma 

■■ ^ " 72 ^ J + T2tarL2(j| 


dip 


(%-) (^) , 



,u 


= 


(^y 


jt 

/ 


sin, (p taxL^'y 


dq> 



-l)2d(p 

(23b) 

(23c) 

d<p 

(23d) 

(23e) 

(23f) 

(23g) 


sets 
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CM 




^2{k)k^dk 


• = 3t2 

J 


ffi(k)k^dk 

'0 



A cos'^ 


(23h) 


The subscript 1 desigojates a longitudinal ccsiiponeiit} subscripts 2 
emd 3 designate the lateral components. Of the remaining subscripts, 

B denotes a shear- flow component; p denotes a potential-flcw compo- 
nent; c refers to the ccaribustible mixture; and b refers to the burned 
gas. Equations (25a) to (23g) have been integrated numerically using 
Simpson’s zule with the following increments for <p: 2° intervals frcan 

0® to 20®j, 5® intervals from 20° to 70°, 2° intervals from 70° to 90°, 
and so forth. Numerical results are listed in table I. 


I 


DISCUSSION OF RESULTS 

Mean-square fluctuation quantities generated by the linearized in- 
teraction of a constant-pressure combustion front with a weak isotropic 
turbulence field satisfying the continuity equation for incompressible 
flow are given by equations (23) in terms of the incident- turbulence- 
velocity caraponents. E^zations (23a), (23b), and (23e) apply only at 
the flame front where the attenuation factors are unity. This restric- 
tion does not apply to equations (23c) and (23d) in the assumed absence 
of turbulent decay processes. 


Velocity Fluctuations 

Potential- flow fluctuations . - The root-mean- square velocities 
associated with the attenuating pressure fields generated in the com- 
bustible mixture (eq. (23b)) and in the burned gas (eq. (23e)) are 
plotted in figure 2. Since these ratios apply only at the flsuae front 
where the ejgjonential attenuation factors are xmlty, they represent 
maxlTmi'm values. As is to be ej^iected from the boundary- condition re- 
quirement of equal pressure fluctuations with differing densities on 
each side of' the flame front, potential-flow velocity ccmponents in 
the burned ges exceed those in the combustible mixture. Both ratios 
Increase with increasing flame- temperature ratio, reaching asyn^ptotic 

values of in the burned gas and unity in the combustible mixture. 

Hot-wire instrumentation will respond to these fluctuation velocities 
as well as to the shear-flow' fluctuation velocities. In view of their 
exponential attenuation characteristics, however, such contributions 
would not be of Importance unless measurements were made at stations 
very close to the flame, that is, within a distance of the order of 
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incident- turbulence scale. For the low flame Mach numbers encountered 
in combustion, the contribution to the hot-wire signal voltage of the 
unattenuated sound waves described in the section FLAME - TURBULENCE 
INTERACTION PROCESS should be quite small. 

Flame- generated turbulence . - The shear flow in the burned gas 
(eqs. (23c) and (23d)) constitutes the flame- gene rated turbulence oc- 
casioned by the presence of approach- flow turbulence. These velocities, 
referred to the incident tiirbulence velocities, are also plotted in fig- 
ure 2. A sli^t ajTq)lification of the longitudinal component occiirs for 
values of t in excess of 7. In the limit, as z becomes very large 
the longitudinal and lateral velocity ratios approach asymptotic values 

of aJzJZ and i *JZ, respectively. 

2 

The diagrams of figure 1 indicate that a pressure wave interacting 
with the flame front csji also bring about a shear flow in the burned 
gas. Although the reflection and the consequent impingement of the 
pressure fields described by equation (23a) upon the flame front are . ... 
possible, any additions to the flame- generated turbulence level throu^ 
the reflection process would probably be negligible because of the at- 
tenuating nature of the pressure field. Thus, contrary to the predic- 
tions of references 6 and 7 that the flame- generated turbulence inten- 
sity should be many times greater than the intensity of the incident 
field, the present analysis indicates that the two intensities are 
about equal. 

It is interesting to note that a stream contraction (ref. 17) in- 
creases the downstream velocity of the mean flow (as does the flame 
front also), while exercising a different selective effect upon an in- 
cident Isotropic turbulence field. For example, with a sevenfold in- 
crease in the downstream velocity of the mean flow, the longitudinal 
velocity ratio (in the absence of decay effects) is 1.01 for the flame 
front and 0.31 for the contraction. The corresponding lateral velocity 
ratios are 0.87 and 2.29, respectively. 


Turbulent Flame Speed 

The hl^er mass-flow birrning rate of a turbulent flame as coaipared 
with that of the corresponding laminar flame is generally described in 
terms of a turbulent flame speed The flame- speed ratio U^j/U is 

generally ass\imed to be equivalent to the ratio of turbulent- to-lamlnar 
flame surface area. 

Calculation of the tiirbulent flame speed reqiilres consideration of 
second-order terms. The local instantaneous normal propagation velocity 
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U + 6U of the distorted flame front into the c embus tible mixture at 
rest is 


CM 

iP- 


CO 

CM 


M 

0 

M 

1 

& 


U + 8U = 


% + 5-t - ^i,c + ^X2 .c^2,c + €x3^3,c 


*JT71 


2 fC 

1 2 j. *r ’ 2 

^2,c 


(24) 


Let Up = U + + U/5\ + . . . where 


0 nhe flan 


contribution to 'the flame speed of order 
ties, for exan^ile, are written as 




represents a steady- state 
The perturbation quanti- 


. Q 

?t = «t 



+ • • • 


uj = + . . . 

1,C 1,C 1>C 

(S) 

The sigjerscript on a fluctuation qiiantity indicates the order of 

the perturbation. As before, SU/U is taken as 


U 


D 


(r - i)“ 2-2 





Substitution of this expression into equation (24) and performing the 
indicated esqjansions yield the foDJLowlng relation, which is correct 
through second- order terms: 



first- second- order 

order terms 

terms 



first- second-order terms 

order 

terms 



28 


MCA ra 3407 


Averaging this equation, yields 



For the limiting case of constant-pressure combustion, F = 0. Thus, the 
ratio of turbulent- to- laminar flame speeds can be written for this case 
as 



The first two terms on the ri^t side of equation (25) represent 
the ratio of averaged turbulent- to- laminar flame surface area. The 
third term is a correlation coefficient representing the transverse- 
velocity- fluctuation contribution to the turbulent flame speed caused 
by flame-front distortion. It is interesting to note that only the 
transverse velocity fluctiiations appear ej^licitly. The second term 
has already been determined (eq. (23g)). The third term may be ob- 
tained from equations (C8) as 


(t - l)tBnS + + (^ + 1)] + i^dq> (26) 

The imaginary term has not been written out inasmuch as it does not con- 
tribute to the integreO.. The results of the indicated integrations are 
listed in table I. 



Equation (25) for the flame- speed ratio Uqt/U may also be written 
in the form 


""t 

= 1 + 


JjLg A .? 

ir 


(27) 


where the flame- speed parameter S is obtained from the values listed 
in table I for equations (23g) and (26). The variation of this param- 
eter with the flame-temperature ratio % shown in figure 3 suggests, 
on the assungption that the flame- frcmt slopes eind ^ gov- 

^2,c ^,c 

em contributions to the right side of equation (25), that the flame 
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front with the hi gh er heat release Is distorted less by a given inten- 
sity of turbulence. For the degenerate case t = 1 (no heat release), 
the parameter S becomes infinite a condition ccanpatible with tfals 
viewpoint. 

The present anstlysis requires that the flame-front slopes, as well 
as the other perturbation quantities, be small. The preceding discus- 
sion suggests that the incremental flame- ^eed ratio (b^ - U)/U may be 
of second order as a result of this restriction to smal .1 flame-front 
slopes. 


Combustion Hoise 


The root-mean- square pressure- fluctuation coefficient 

which applies directly at the flame front, is plotted in figure 4 for 
the limiting case of constant-pressure combustion. The pressure fluc- 
tuations are a measure of the random noise generated by the interaction 
of the flame front with the incident turbulence. In acoustical measure- 
ments the noise level in decibels is usually given with respect to a 
reference pressure of 0.0002 dyne per square centimeter (ref. 18), which 
corresponds approximately to the pressure angplitude of a plane sound 
wave of minimum, audible intensity at a frequency of 1000 cps. The noise 
pressure level in decibels is defined by the relation 



74 + 20 log]_Q 

(28) 

where the pressiire fluctixations in the combustible mixture are given in 
dynes per square centimeter. 

Equations (23a) and (28) indicate that the noise level should be 
particularly dependent upon flame speed. Propane-air and acetylene- 
air combustion, which are characterized by a low flame ^eed etnd a 
fairly hi^ flame speed, respectively, will be considered for illus- 
tratlve purposes, Pertinent data for these flames at maxlmum-fl ame- 
speed and stoichiometrie conditions for an ambient tenperature of 25° C 
and a pressure of 760 millimeters of mercury are given in table H. The 
adiabatic equilibrium flame tenperatures, at which the total enthalpy of 
the fuel and oxidant equals the total enthalpy of the products of reac- 
tion, were calculated using the procedvire of reference 19. (Total en- 
thalpy includes the chemical contributions to the internal energy.) 
Flame-speed data were obtained from references 20 and 21. 



V^c 

noise pressxire level = 20 logio yo odo2 
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If the flame-front turbulence intensity 
equal to 10 percent, noise pressure levels of 59 and' 81 decibels are 
obtained for propane-air flames and acetylene-air flames, respectively, 
under conditions for maxiinnm laminar flame speed. At an approach- flow 
velocity of 1225 centimeters per second, which is in the range of ve- 
locities usually encountered in combustion experiments, the correspond- 
ing Intensity of the approach-flow turbulence would be about 0.3 percent 
for the propane-air mixttire and about 1 percent for the acetylene-air 
mixture . 

Thus, the pressure fluctuations generated at the flame front when 
the incident turbulence is of low intensity, although small as compared 
with ambient pressure, are apparently of fairly hi^ acoustical inten- 
sity for constant-pressure combustion. Because of the exponential at- 
tenuation of these pressure fields, the "far-field" acoustic intensity 
(the intensity at distances very far from the flame front) approaches 
zero. For cases other than constant-pressure ccmbustion, a finite 
"far-field" Intensity is obtained. 


CONCLUDING REMAEKS 

The present linearized analysis has treated the interaction of a 
field of isotropic turbulence with a free flame front under constant- 
pressure combustion conditions with no turbulence decay processes or 
heat- release fluctuations. The interaction produces an anisotropic 
turbiilence field in the burned gas which has axisymmetry about the 
main-stream direction. Contrary to the results predicted by several 
current theories of turbulent flame speed, the flame-generated turbu- 
lence velocities caused by approach-flow turbulence do not differ 
greatly from the turbxilence velocities of the incident field. 

The incremental flame-speed ratio (u^ - U)/U as obtained from the 
present analysis is a second-order quantity consisting of two parts: 

One part represents the root-mean- square area of the turbulent flame 
fronts the other represents the contribution of the transverse veloc- 
ity fluctuations which result from the flame-front distortion. The 
flame-speed ratio ^/U for a given level of Incident- turbulence in- 
tensity U is found to decrease with increasing heat-release 

rates (increasing values of r). 

Random pressure fluctuq,tions generated in both the combustible mix- 
ture and the burned gas, although small compared with ambient pressure, 
give rise to appreciable noise levels (59 to 81 decibels) directly at 
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the flame front even for very moderate intensities of approach- flow tur- 
bulence (flame- front turbulence intensities of 10 percent). For the 
limiting case of constant-pressure combustion, the pressure waves at- 
tenuate exponentially with distance from the flame front, so that the 
"far-fleld" intensity approaches zero. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, January 25, 1955 
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APPENDIX A 


SYMBOLS 

The following synibolB are used in this report: 

A magnitude of two-dimensional vortlclty-wave velocity 

vector in combiistible mixture 

A two-dimensional vorti city-wave velocity vector in com- 

bTistible mlxtiire 


a 

b 

C 


"P 

°v 

Dl 

d 

El 

Fij(^,t) 

f 

G 


H 


speed of sound 

coefficients defined in eqs. (4) 

constant defined in eqs. (7) and (9) 

amplitude of combustible mixtiure shear-wave ccai 5 >onent 
parallel to x^^^-axis 

constant defined in eqs. (7) and (9) 

specific heat at constant pressure 

specific heat at constant volume 

U/(U - V), eqs. (4) 

constant defined in eqs. (7) and (9) 

-V/U, eqs. (4) 

spectral tensor function 

constant defined in eqs. (7) and (9) 

(G^[_ + 1 G 2 )/A sin Q> 

an^Jlitude coefficients of shear-wave longitudinal com- 
ponent in burned gas 

(H^ + 1 H 2 )/A sin q> 

an 5 >litude coefficients of flame-front displacement 
hio groupings defined in eqs. (B16) 
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in 

I 

•t) 

o 






an^illtude coefficients of shear-wave transverse ccmi- 
ponent in humed gas 

amplitude coefficients of pressure wave in burned gas 


j(l) j(2) j(l) = — ^ Jn; s — ^ Jsf eq. (lOd) 

Pc^ Pc^ 

Ki,Kg,K 3 ,K^ coefficients defined in eqs. (4) 

k magnitude of wave-number vector k 

k wave-number vector of shear wave in ccmbustlble mix- 

ture with components k', k' in x ,y- coordinate 

12 c 

system, with coniponents k-j_i k£, kg in j,, 

*2,c^ X 3 ^g- coordinate system 

L^_,L 2 amplitude coefficients of density associated with 

"shear entropy wave in burned gas 

M flame-front Mach number, M s U/a.^ 

m^ a^/yU^, eq. (B16) 

N (% + iH 2 )/A sin <p 

Nj_^N 2 an 5 )lltude coefficients of longitudinal velocity com- 

ponent associated with pressure wave in burned gas 

n exponent used in representation of laminar flame speed 

as function of cambustlble-mlxture static temperature 

P^,P 2 anplltude coefficients of lateral velocity con^onent 

associated with pressure wave in burned gas 


P 

P' 

Q 

E 




static pressure 
statlc-pressTire perturbation 

term not contributing to the integral in eq. (26) 

(e^^^ + R^^^)/a sin <f> 
gas constant 

amplitude coefficients of pressure wave in combustible 
mixture 
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r(1),r(2) 

r 


S 

T 

Tij(r,t) 

To 


= R 3 _/')fM?, R^^^ = Rg/yM^^ e<l- (lOc) 
separation vector 

constants used in representation of lamlneo* flame 
^eed as function of cambuBtible-mlxture static 
tenperature 

flame-speed parameter, eq. (27) 
static ten^erature 

velocity correlation tensor for homogeneous turbulence 

stagnation teurperature 

time 


U 


laminar or fundamental flame speed 


Uj 


mean turbulent flame speed 


u' longitudinal conponent of velocity perturbation 

'^l c^^2 c velocity perturbation con^jonents in conibustlble mlx- 

^ ^ ’ ture parallel to Xg^^; X3 ^p coordinate axes, 

respectively 


^l,b^^2,b^^3,b 

V 


velocity perturbation components in burned gas parallel 
^l,b'’ ^2,b-' *3,b coordinate axes, respectively 

mean velocity of burned gas 


'T 


mean velocity of biimed gas in turbulent combustion 


V' 


lateral component of velocity perturbation 


W 


(W]_ + IWg) /a sin <p 


amplitude coefficients of longitudinal velocity com- 
ponent associated with, pressure wave in combustible 
mixture 


w' lateral ccmiponent of velocity perturbation (component 

parallel to plane of unperturbed flame front and 
normal to u' and v' components) 




amplitude coefficients of lateral velocity component 
associated with, pressure wave in combustible mixture 


A 
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Xc 




^2,b 


^3,b 

*l,c 


X 


2,c 


*3,c 

y 




coordinate in x^,y-system measured In direction of 
unperturbed flame-front travel relative to which 
burned gas is at rest 

coordinate in Xj,,y-system measured in direction of 
vmpertiirbed flame-front travel relative to which 
combustible mixttire is at rest 



rection of unpertiirbed flame-front travel relative 
to which burned gas is at rest 


coordinate orthogonal to Xjj_ ^ and Xj and making 
angle 0 with y- coordinate 

coordinate orthogonal to x^_ ^ and Xg ^ 

coordinate in x^ c^*2 c^^3 measured in di- 

rection of unperturbed flame-front travel relative 
to which combustible mixture is at rest 

coordinate orthogonal to X]_ ^ and X3 ^ and making 
angle 0 with y- coordinate 

coordinate orthogonal to and 

coordinate orthogonal to x^. and x^ 

random amplitude vector of shear- field Fourier 
component 


z 


r 

r 

A 

« &0 


coordinate orthogonal to x^, and y 
grouping defined in eqs. (B18) 


1 

(r - l)nr 2 Tc^” 2 m 



ratio of specific heats 

- l)^tan'^cp,+ 2 t(t^ + 2t - l)tan^:p + (t + 1)^ 
grouping defined in eqs. (B24) 
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grouping defined in eqs. (B2l) 

varialQle upon which pressure wave depends, eqs. (?) 
variable upon which pressure wave depends, eqs. (7) 


® 1»®2 • • 


0 


A 


V 





p' 


a 


T 

$ij(k,t) 

q> 


s 


‘10 


angle between polarization plane of incident shear 
wave and XQ^^Q,X 2 ^£j-plane 

coefficient defined in eq. (3e) 

flame- front displacement 

flame-front displacement velocity 

flame-front slope with respect to X 2 ^j,- coordinate 

flame-front slope with respect to X 3 ^j.- coordinate 

flame-front slope with respect to y-coordinate 

static density 

static- density perturbation 

k£ut + kgy, eq. (lOj) 

flame-ten 5 )erature ratio, Tg^-| 3 /T 0 ^j, 

spectral density tensor 

angle between wave-momber vector of incident shear 
wave and direction of unperturbed flame -front 
travel, tan = iSg/kj 



scalar amplitude function defining spectral density 
tensor 


A 


Subscripts; 

1,2,3 


orthogonal coordinate designation 
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b 

c 

cr 


P 

s 

Superscripts: 

(S) 

* 


burned gas 
ccmbustible mixture 
critical 

designates order of steady-flow quantity 
potential-flow velocity ccanponent 
shear- flow velocity component 

designates order of fluctuation quantity 
denotes complex conjugate 

denotes fluctuation quantity except where otherwise 
specified 


t 



38 


NACA TH 3407 


APPENDIX B 


ATTHSUATING-WAVE SOLUTION FOR SINGLE-WAVE INTERACTION 
set of equations; 


+ y) - ®2Gi = - ®2 ^ sin «!> - BgWi + B3R1 

+ ^] + B1L2 - B2N2 - B2G2 = - B2W2 + B3R2 

Ji (i' r) " ^ 

Jg^l. - L2 - K1N2 - K1G2 = - %V 2 + K2R2 - (%- ^4)^2 

i J^ + Li + (Di . l)Hi - DiNi - DiGi = - Wi - A Bin (P + ^ 

i J2 + Lg + (Di - 1 )H 2 - B1N2 - D1G2 = - W2 + i R 2 j 


>(B3) 


Pn + Ii = Xi - ^ A Bin q> + p- EiHi 


^2 




(B4) 


tr ^ ) % " '^b®2 = m^(-i‘b'^l + ^b^2) 
<ib% - - ^bJ2) 


(B5) 
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^bPl - fb 


- <^^^2 = “^"'<>bJ2 

"^l(u ? v)®2 = ^2^: 

-^l(u^)^ = ^2ll 

fc^i - d^W2 = (f^Ri + b^Rg) 

yM 

i-^A + 


(B6) 


(B7) 


(B8) 


^c^l ~ ^c^2 ~ ^c^2 

dcXi + f^Xg = - C^lj 

H2_ - - A sin (p = AR^^^^ 

^2 “ ^2 = 

From eqxiations (B5), (B6), (B8), (B9), and (BIO), 

% = 

No = 


^1 ~ ^l(^2'^l ■ 

^2 = ^l(^3'^l ■'■ ^2*^2 )/ 

i " (-114^1 - 1^J2)\ 

'2 = - ^ 4 ^ 2 ^ J 

+ hgRg) 

+ hyRg)/ 


^’l = 

?2 = 

^2 = ^(“^ 8^1 
Xi = h0(-hgR2 + 1110^2 

^2 " ^(-^10% " 119^2)/ 


(B9) 


(BIO) 


(Bid.) 


(B12) 


(B13) 


(B14) 
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= AR-j^ + A sin + hg 

Hg = ARg ■*" ^(“^8^1 ■*" ^7^2^ J 


where 


^1 = 


f |/Df + df ^ ^ rM2(f2 + d2) 




1^2= le - 


h3- 




Dn + ^b^b 


h4= c^d^ 


^5 = 


Dn 


h^= - b d 

7 c c c 


hQ= b^f^ + d^f^ ^ 
^9 ~ '^c'^c 


1^10- 


2«. 
m s 


%/ru^ 




From equations (B4) , (B7), (B12), (B14), and (B15), 

G. « 


where 


\ « a^Rl - cCgRg - a3J^_ - + a^A sin (p^ 

Gg = ^ 2^1 "* ^^^^2 } 

%. 

^2 l-^2\ 1 

= k]5^ ^ 

hp r hpEn I 

° ^ 


^2 






“4 = lull's 


k’ ^ 

“ 5 * 

^1 


J 


(B15) 


(B16) 


(B17) 


(B18) 


kb' 
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Then, from equatlcjn (B2), 

r 1 

Li = Jj[l - ^ - K-^Chihg - a.3)J + “45] ^ 

Rl[Ki(h6h7 - ctq^) - Kg + (Ki - K4)AJ + 

■*■ °'2 ^J <p[Ki (1 “ “5)3 “ 

Lg = - J3_[K2^(h2h3 + a^)] + J2|l - •i - Kj^Ch^hg - ct3)J - R3^[K2^(h0h0 + ag)] + 

Rg[Ki(h0h7 - ai) - Kg + (Ki - K4)A] J 

(B19) 

The various disturbance amplitude coefficients of equations (lO) 
have now been obtained in teims of the coefficients Rj_, R2# J2 
and the parameters A sin <p and t*. From, equations (Bl) and (B3), 


JlSi + Jg6g + RjB3 + R264 = £ 5'^ 


-Jj^fig + ^2®! “ ^ 1®4 ^E®3 ~ ^ 

JlCg + J2®7 ■*■ %®9 = ®10 


> 


-J167 + JgBg - RjL6g + RgSg = Oj 


(B20) 


where 


®1 - 1 + - (^1^1 + Bg)(hihg - a3) 

®2 ~ ^®1^1 ■*■ ®2^(^1^3 “4^ 

®3 = (Bi^i + B2)(hgh7 - aQ_) + B2 _(Kj_ - K 4)A - B3 - B^Kg 


®4 ^ (®1^1 ■*■ ®2HB6b.0 + 02) 

6g * B2K3T * — A sin [([BjKi ®2 ) ( ^ “ “5 3 

eg s 1 - + Ki)(hih2 - 05) 


>(B21) 


67 2 + Kj_)(h3^h3 + a4) 

Eg S (D^ + Kj) (hghy - o^) + (D^ - 1 )A - ^Kg + + (Ki - K4> 

Bg S (D-l + Kj^)(hghg +0,2) 

®10 ~ ^3'*'' “ ■*■ + K3^)(1 - ag) y 
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or 


where 


= 


^2 = 


*2®6 “ ®1®7 

~ *4^7^ ~ ^2(^2^8 ~ 
®2®6 " ®1®7 


Ri = 


^2^6 ” 


'3°5 


®2®4 “ 

6i8g - 8364 
2 6264 - 8 t8 , 


'1^5^ 


2 «i(®i®9 - ®4®6^ *3^®2®6 “ ®i®7^ " ®2^®4®7 " ®2®9^"^ 

&2 2 ®i(«ie0 - CjCg) - 64(6266 - ej_67) - 63(636^ - SgCg) 

^3 - ®5(®2®6 “ ®1®7^ 

84 - ®e(®l®9 “ ®4®6^ + ®a(®2®6 - ®1®7^ ” ®?(®4®7 " ®2®g) 

®5 ® '6<V8 - '3'e) - '7<‘3'7 ' ' 's^Vs ' V?) 

^6 ® *10^®2®6 " 


(B22) 


(B23) 


(B24) 


Equations (BU) to (B15), (B17), (B19), (B22), and (B23) provide the 
formal solution of equations (B1) to (BIO). 
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APPENDIX C 


ATTENUATING-WAVE SOLUTION FOR CONSTANT-PRESSURE COMBUSTION 

If terms of order M^ are retained, equations (l) provide the 
following relations for the unperturbed flow quantities: 




= 1 - - i)m2 + . . . 




°P^S,c _ 1 

r\ » 




2 ir - 1)m2 


u2 _ 

a2 ' t[2 - (r - 1)(T - 1)m2 - ir^ - 1)t(t - l)M^ + . . .] 

With these relations, equations (4) take the form 

B2_ = r^M2, Bg = 2 tM 2, Bg = 1 + [l + rC'T - 1 )]m2 

D^ = (t - 1)m 2], El = (t - l)[l + tM^] 


Ki = K4 = (r - i)t^. Kg = ^^^[1 + K - i)m2]. 




> 


Tf - i 4. r - 1 

K3 ^ + 2 

(Cl) 


and equations (9) may be written 


f2 = k'2 


(k^2 _ ^ jj^2 




(C2) 
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If only the leading terms in powers of are retained, equations 

(B16) and (B18) of appendix B provide the following relations; 


hihg - hj^hg = 


T^tan^ X tan 

jr2Vn ^ ^1^3 - ^1^4 = „2 


VM^(l + x‘^tenr<f)' 


hghy = hghio = 


^6^8 ^ ^^9 = 


rM^(l + T^tan^cp) 
sin <p cos (p 


rM^ 


(C3) 


and 


ai = 


T sln^ 


YMt 


— (l - (t - l)tan^<^ 


tto = 


T^sin^y tan q> 
yM^ 

T^tan^cp 

yM^(l + T^tan^fp) 

T^tan^(p 

^ rM^(l + T^tan^q)) 

=? t[i - (t - l)taji^«^ 


(C4) 


From these relations: 


1 . tan q> - l)sln2cp(‘’^ tan^qj - l) 

hih3 + = • -2 ^ ^6^17 - «1 




tiqhg - tt3 = 0, 

From equations (B21) and (C5), 


, , , sin «p COB q>(l + T^tan^q)) 

hghQ + ag 


(C5) 


_ _ sin q> cos ^(l + T^tan^q>) 

®1®9 " ®4®6 = T3 

yM^T 


®1®8 " ®3*6 ~ 


(t - l)sin^q>(T tan^q> - l) 


yM^T 


> 


®3®7 “ ®2®8 = ®2®6 " ®1®7 = " } ®4®7 “ ®2®9 = 

yM^ ^ 


(C6) 
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aad from equations (B2l), (B24), and (C6); 


5]^ = *P + l)(l + T tan^<p) - 2 tan^<p(T - 1 )t(t tan^tp - l)| 

rM^T 

82 = |(t - 1)(t tan^q> - l) + 2 (t + 1)t(t tan^<P + 1^ 

63 = - tan cpjj' - 2 (t - i)(t tan^ip - l)A sin «p] 
e -(t - l)tan <P sin^(T tan^<(> - l) 

c ('T + l)('^ tan^<P + l)sin^<P 
85 = - - (t - 1) (t tan^ - 1)A sin 




J 

(C7) 


Since for constant-pressure combustion = Rb/tM^, 

and so forth., the coefficients r(^), and j(2) 

are obtained 

from eqiaations (B22), (B23) , (C6), and (C7). With these coefficients 
determined, the remainder are obtained from equations (B7), (B11) to 
(B14), (B17), and (B19). The amplitude coefficients are 
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“ X 'C-[t2(t2 + 2t - l)tflji2iy + t(t 2 + 3^)j 

r( 2) .= j(2) „ _ ^ ^ !{-c 2 . i)(-c2tan4^ _ + _ i)tan*p - -ar^taii^q) - tCt + ij] 

“ ■ 2 tan^«P - l)^A sin ^ + t2[t(-i - l)®tan*qi + (t^ + 2t - l)tan^* + aj 

Ll - - f^, L2 - 0 

Hi = ^ j tan^^t - l)[(-c + 1 )(t tan2((i + 1) -^(t - l)tan®(p(T tan^ip - 1 )]a sin f - 

t^^tan®? ^(-u - l)^tan^^ + tC-c + l)tan®^ + t]+ (t + 

Ho = ~ Vl"^ tan^^ - + 1 )(t tan^<p + 1) + (-t - 1 )(t tan^p - 1)1 a sin p + 

A(1 + T^tan^p) (_ 

i2['c( 2^ . i)tan^P - (t - l)tan2p - «} 

On ■ — ^ 5 — ^rcv + 1)2 (t tan^ + 1 )^ + tOv - l)^taiSp(T tan^ - 1 )^JA aln p - 

A(l+t‘taa‘^) I"" -' 

t(T - l)tan^[T(T^ + 2 t - l)tan^p + (t® + 

Qg = . f . i)(t tan^p - 1 )(t tan^p + l)A slnp + - l)tao*P - St^tan^p - t{'c + 


> (C8) 


?1 = - Mg, Pg = 


" T tan tp' ^2 " tan 9 


1 . l)(, tan^p - 1)[^T - 1) Wp + et + 1 )]a Slap - [t^CSt - l)tan®p + t(t + 1)]^ 

^ _ 7 ^ T - l)tan P sin p + t[t tan^P + 


% “ 'fg, ' ^*1 


El = 2 + l)('T tan^p + 1)A sin p+ fT(3r - l)tan®p + (t + 


3493 



3493 


NACA TN 3407 


47 


For constant-pressure combustion = ^ 1 ^ = k^Ut + k^y, 

"Ic = \ = ^2[-^ + '' = ^i^'c ^2^^ 

i|f = — — — + kly, and a = k'Ut + kly- At tlie flame front where x = Ut 

and Xij = (U - V)t, ''Ic ~ ~ attenuation factors e ^ and e ^ 

are unity^ and the arguments of the dlsturhance waves have the form 
(kj^Ut + k^y). 

The proced^lre of reference 23 is utilized in reference 22 to treat 
a similar interaction problem for consideration of flame-front stability. 
It is assumed that there are first-order perturbations in the laminar 
flame speed U (in the present notation). Reference 22 presents results 
only for the special case where the plane of polarization of the inci- 
dent shear wave is in the plane Xj^ ^ - Xg = 0°) and the wave-number 
vector is parallel to the unperturbed flame front (<p = 90°). The re- 
s-ults of the present analysis were compared with those of reference 22 
for the case of an absence of first-order perturbations in the laminar 
flame speed and for the special case of 0=0° and <P = 90° without 
heat- release pertiirbations. Althou^ agreement as to sign and magnitude 
is obtained for the shear-field velocity components, apparently differ- 
ences in sign occur for the potential- field veloclty-coniponent ampli- 
tudes as indicated in the following table: 



Present analysis 

Analysis of 
reference 22 

U 

-A 

-A 

I 

V 

U 

iA 

-iA 

U 

A 

-A 

u 

iA, 

iA 


Consideration of conservation of momentum across the flame front indi- 
cates that Vp^(,/U and Vp -^/U should be of the same sign. The pres- 
ent analysis is in agreement with this consideration. 
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TABLE II. - BIDROCAEBON-AIR FLAME DATA AT SPECIFIED COHDITIONS 


(Pc = 760 mm Hg; T^ = 25° C) 
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(b) Interaction for wave-ntmiber vector inclinations generating noncoalescing 
cylindrical sound waves . 


Figure 1. - Wave formation arising from interaction of flame front with vortlcity wave. 
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Root-mean-aquare velocity- 
fluctuation ratio 
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Figure 2. - Effect of flame-temperature ratio on shear-flov and potential- flow velocity 
fluctuations . 
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Flame-temperat\ire ratio, t 

Figure 3. - Effect of flame-temperature ratio on flame speed for incident 
Isotropic turt)ulence. 



Root-mean- square noise pressure-fluctuation 
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Flame-temperature ratio, T 


Figure 4. - Effect of flame- temperat\are ratio on random 
pressure fluctuations generated at flame front. 
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